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ABSTRACT 

The uncertainty in the determination of the momentum and scalar fluxes remains 
one of the main obstacles to accurate numerical forecasts in low to moderate wind 
conditions. For example, latent heat fluxes computed from data using direct covariance 
and bulk aerodynamic methods show that there is good agreement in unstable conditions 
when the latent heat flux values are generally positive. However, the agreement is 
relatively poor in stable conditions, particularly when the moisture flux is directed 
downward. If the direct covariance measurements are indeed accurate, then they clearly 
indicate that the bulk aerodynamic formula overestimate the downward moisture flux in 
stable conditions. As a result, comparisons of the Dalton number for unstable and stable 
conditions indicate a marked difference in value between the two stability regimes. 

Investigations done for this thesis used data taken primarily at the Air-Sea 
Interaction Tower (ASIT) during the Coupled Boundary Layers and Air-Sea Transfer 
(CBLAST) Experiment 2003 from the 20-27 August 2003. Other data from the shore 
based Martha's Vineyard Coastal Observatory (MVCO) and moored buoys in the vicinity 
of the ASIT were also incorporated. During this eight day period, the boundary layer was 
often characterized by light winds, a stably stratified surface layer and a swell dominated 
wave field. Additionally, the advection of warm moist air over cooler water resulted in 
fog formation and a downward flux of moisture on at least three occasions. Therefore, a 
primary objective of this thesis is to present a case study to investigate the cause of this 
shortcoming in the bulk formula under these conditions by examining the physical 
processes that are unique to these boundary layers. Particular attention will be paid to the 
behavior of the Dalton number in a stable marine atmospheric boundary layer under foggy 
conditions using insights derived from the study of fog formation and current flux 
parameterization methods. 

Thesis Supervisor: Dr. James B. Edson 
Title: Senior Scientist 
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1. INTRODUCTION 


1.1 Background 

Over the past several years, WHOI researchers and their colleagues have made 
significant progress in obtaining measurements of turbulent fluxes in the marine 
atmospheric boundary layer (MABL) through recent field experiments such as TOGA 
COARE (Fairall, Bradley et al. 1996), RASEX (Mahrt, Vickers et al. 1996), MBL 
(Grachev, Fairall et al. 2003), COOP (Hara, Bock et al. 1998), and FASTEX (Hare, 
Persson et al. 1999). The majority of these measurements were collected in near-neutral 
to slightly unstable conditions. The culmination of these efforts resulted in the state of 
the art TOGA COARE 3.0 (TC 3.0) bulk aerodynamic (BA) parameterizations (Fairall, 
Bradley et al. 2003). This BA parameterization has been widely accepted and proven as 
the field standard for turbulent flux measurements in the MABL, particularly for unstable 
conditions. The current parameterization for stable conditions is based on the SHEBA 
(Persson, Fairall et al. 2002) experiment which was conducted over ice. However, 
several of the aforementioned experiments collected data in stable conditions when the 
direct covariance (DC) and TC 3.0 latent heat flux (LHF) differ significantly. 

Specifically, this data indicated that when stable conditions existed, particularly in fog, a 
downward (negative) flux of moisture to the sea surface was frequently observed. A 
negative surface flux means that the ocean is gaining heat and the atmosphere is losing 
heat. As can be seen in figure 1.1, TC 3.0 tends to overestimate the LHF when compared 
to DC measurements in these conditions. Additionally, with the exception of RASEX 
and MBL, these experiments did not measure flux profiles and neither RASEX nor MBL 
was set up to examine mean and flux profile relations for scalars such as moisture and 
heat. As a result, we still rely heavily on overland measurements such as those made 
during the Kansas (Izumi 1971) and Minnesota (Champagne, Friehe et al. 1977) field 
campaigns for flux profile relationships, in addition to the above mentioned relationships 
over ice, for unstable and stable MABLs. 

Recently, flux profile relationships in the MABL were measured during the 
GasEx and FAIRS experiments (Edson, Zappa et al. 2004; Hintsa, Dacey et al. 2004). 
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The data from these experiments exhibited good agreement with overland experiments 
for unstable conditions, but were unable to evaluate results in stable conditions due to a 
lack of sufficient data. In addition, other complications arose that are inherent to any 
experiment where measurements are taken strictly from ship or buoy-based 
instrumentation platforms. These issues arise from the requirement that corrections must 
be made to remove the effects of platform motion, localized 

Figure 1.1 CBLAST Latent Heat Fluxes August 20-27, 2003 



heating, flow distortion (particularly the ship-based GasEx measurements) and instrument 
contamination due to the adverse effects of the marine enviromnent such as sea spray, 
corrosion and wave action. 

The Office of Naval Research (ONR) funded the Coupled Boundary Layers and 
Air-Sea Transfer Low Wind (CBLAST-Low) experiment to further the progress already 
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achieved in the previous field experiments under light to moderate wind conditions. The 
experiment was designed to measure fluxes and their associated profiles so that every 
term within the turbulent kinetic energy and scalar variance (SV) budgets were calculated 
as described in section 2.3. An optimal suite of instrumentation deployed from aircraft, 
ships, buoys, and fixed-towers acquired the most complete and comprehensive data set 
ever collected in the MABL. Simultaneous meteorological and oceanographic 
measurements provided scientists with a unique opportunity to fully investigate the 
complex physical processes which drive the energy exchange between the ocean and 
atmosphere. 

One of the primary objectives of the knowledge gained from the CBLAST 
experiment is to aide in the development of improved parameterizations for possible use 
in the Navy's mesoscale model, the Coupled Ocean-Atmosphere Mesoscale Prediction 
System (CO AMPS). Despite its name, CO AMPS is not, as of yet, truly coupled. This is 
mostly attributable to the challenges alluded to previously. In particular, as with most 
models, accurate boundary conditions are essential to proper initialization and subsequent 
proper operation of the model. The complexities of the marine environment make this a 
formidable task and are introduced in further detail here. 

1.2 Challenges 

The differences between the land and marine environments are quite distinct. The 
large heat capacity of the ocean, the relative importance of advective processes, and the 
complexity of the air-sea interface make accurate forecasts of sea surface temperature 
(SST) challenging. The presence of waves introduces complications that invalidate the 
application of parameterizations derived over land to grid points near the surface. For 
example, modification of the near surface turbulence by the waves impacts the fluxes 
driven by the turbulence. This invalidates one of the primary assumptions of commonly 
applied scaling laws used in the parameterization of these fluxes known as Monin- 
Obukhov similarity described in section 3.1. In addition, the effect of fog in a stable 
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boundary layer on the turbulent flux profiles over the ocean is, at present, relatively 
unknown. With these complications and despite the progress made in recent years, the 
operation of truly coupled mesoscale prediction models remains a challenge. 

1.2.1 Sea Surface Temperature 

In the case of CO AMPS and other models used over the ocean, an accurate 
prediction of the evolution of SST is a paramount factor. Accurate sea surface 
temperatures are required to yield accurate heat flux estimations from bulk 
parameterizations. Likewise, accurate flux estimates are necessary to provide realistic 
sea surface temperatures and to keep the predictions from "running away", particularly in 
regions or times when data assimilation is an issue. This proves quite difficult as SST 
may have rich time and space variability in low winds with large spatial gradients from 
onshore to offshore. Additionally, the oceanic boundary layer (OBL) changes in 
thickness and temperature from onshore to offshore, depending on, e.g., whether 
upwelling or downwelling favorable winds are present and whether the OBL is 
influenced by the bottom boundary layer (Weller 2004). Due to the huge heat capacity of 
the ocean as compared to the atmosphere, latent and sensible heat fluxes are strongly 
correlated with the slowly evolving SST field (Farrar 2004). If there is a positive LHF, 
the effect of evaporation of the sea surface creates a cool skin effect. If the LHF is 
negative, the effect of condensation on the sea surface creates a warm skin effect. Either 
process produces a temperature gradient in the uppermost millimeter of the sea surface of 
a few tenths of a degree Kelvin. Bulk measurements of sea temperature cannot account 
for this effect so other methods must be employed such as infrared radiometry. However, 
if fog is present, these may give inaccurate measurements. 

SST can also be affected by surface films or surfactants composed primarily of 
phytoplankton. Bands of these surfactant slicks are related to SST anomolies. These 
slicks reduce the surface tension of the sea surface thereby modulating physical transfer 
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processes. These slicks erode with increasing winds and so tend to be more prevalent in 
low wind conditions such as those experienced during CBLAST (Frew 2004). 

1.2.2 Wave Effects 

Waves have a direct effect on the momentum and turbulent kinetic energy (TKE) 
fluxes. The physical argument is that momentum can be carried by both tangential stress 
(i.e. shear stress) and normal stress (i.e. form drag) created by the interaction between the 
pressure and wave fields, which causes interaction between the pressure and velocity 
field away from the surface. The interaction drives fluxes and flux divergences that 
appear as source/sink terms in the TKE and momentum budgets (Janssen 1999). 

However, the equivalent of a normal stress does not exist in the heat and moisture budget 
equations due to the absence of pressure terms. This translates to the absence of a wave 
induced heat or moisture flux over waves (Edson, Zappa et al. 2002; Stull 1988). Instead, 
the effect of waves on heat and moisture fluxes is an indirect result of the waves 
modifying the turbulence responsible for transporting these passive scalars. As a result, it 
is assumed that the scalar fluxes are less affected by waves. 

The form drag of the longer, fastest moving waves can actually impart momentum 
back into the atmosphere (Grachev, Fairall et al. 2003). This effect of form drag on 
momentum exchange and the near surface velocity profile is a function of wave age and 
stability. In light winds, old seas, and stable conditions, the wave generated winds can 
produce a low level jet at 50-100 meters in the MABL (Sullivan 2004). Therefore, 
surface wave effects are essential for coupling the MABL to the OBL in any model 
(Edson 2004; Hristov 2004; McWilliams 2004; Sullivan 2004; Vickers 2004). This 
region where the wave effects have the most pronounced affect on the atmosphere is 
called the wave boundary layer (WBL). At present, the WBL is, in quantitative terms, 
poorly understood. As a result, no generally accepted definition of the height of the 
WBL exists. Some modelers assume that it is limited to the region « 1, where <J H is 

the significant wave height of the dominant waves and z is some reference height. Some 
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recent field campaigns have shown that some terms like the pressure transport term in the 
TKE budget equation are influenced by waves up to heights where k p z~ 2 where k p is 

the peak wave number of the dominant waves, suggesting a much thicker WBL (Edson, 
Zappa et al. 2002). Further knowledge on the WBL gained through the CBLAST data set 
should provide valuable insight into the details of the physics related to this region. 

1.2.3 The Stable Boundary Layer 

In the unstable MABL, current parameterization methods perform relatively well 
except for very near the surface in the WBL. Like its ocean counterpart, the greatest 
challenge in micrometeorology is to improve our understanding of the stable MABL, 
which is often characterized by light winds and fog. Shallow stable boundary layers are 
common in coastal regions and are often characterized by a turbulent flow that is very 
weak or collapsed entirely (Sullivan 2004; Vickers 2004; Mahrt 2004). In fact, the 
boundary layer is often indefinable in very stable conditions. Additionally, the low winds 
mean that there will be less wind generated waves, which allows swell to dominate the 
wave field. This was a common feature during CBLAST. As previously mentioned, 
these swell waves, if they are moving fast enough, can impart positive momentum flux to 
the boundary layer (Grachev, Fairall et al. 2003). 

1.2.4 Fog 

The requirement to accurately predict fog and low-level clouds are big problems 
for all mesoscale models (Wang 2004). The formation process of fog is different from 
that of a cloud. In cloud formation, the air is saturated by adiabatic cooling due to a 
falling pressure in rising air parcels. Fog occurs in the lower atmosphere within a few 
meters, few tens of meters, or at most hundreds of meters of the surface, and is 
commonly driven by horizontal advection and vertical mixing of temperature and water 
vapor (Binhua 1985). This is especially true in the marine environment where advective 
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processes play a key role. However, whether or not the exchange of heat and mass leads 
to fog formation depends on the synoptic state of the atmosphere and the underlying SST. 
Clearly, air-sea interaction plays a key role in fog formation over the ocean. As a result, 
fog presents a multitude of difficulties to scientists and modelers due to the complex 
thermodynamic processes that lend to its formation. 

How fog affects the momentum, heat and moisture flux profiles is largely 
unknown as few if any field experiments have emphasized turbulence measurements in 
fog (Welch and Wielicki 1986). What observations have shown (mostly over land) is 
that fogs are not uniform, but have high frequency variations in most measured variables. 
They also have quasi-periodic oscillations in measured properties such as temperature, 
liquid water, visibility, wind speed, radiation, and turbulence generation. Periods of 
oscillation range from 5-30 minutes depending upon the cause of the oscillations. Some 
attribute the quasi-periodic oscillations to gravity waves propagating at the top of the 
boundary layer. Longer period oscillations may be caused by a fluctuating balance 
between radiational cooling at the fog top and turbulence generation. Shorter period 
oscillations may be caused by advection of fog cells past the observation site. If wind 
direction and wind speed vary slowly with height, then two dimensional cells are 
typically formed with longitudinal bands elongated along the wind direction (Welch and 
Wielicki 1986). 

There are three primary processes which influence fog formation: cooling, 
moistening, and the vertical mixing of air parcels with different humidities and 
temperatures. In the marine environment, warm, moist air moving over cooler water 
results in sea fog, (sometimes called advection fog), while cool, dry air moving over 
warmer water generates steam fog (sometimes called radiation fog). Historically fog has 
been classified by these two main categories: radiation and advection (Lundquist 2000). 
The dominant fog species present during the case study period was advection fog. 
Therefore the remainder of this discussion will focus only on advection fog caused by the 
advection of warm, moist air over the relatively cool water south of Martha's Vineyard. 
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One school of thought is that the cool water in contact with the air causes the air 
to cool to its dew point and condense, forming fog. However, this layer of air by contact 
is generally very thin, and condensation, if any, is thought to occur in a very thin layer, so 
that the fog formed may also be shallow (Binhua 1985). However, advection fog often 
exists in thick layers. Therefore, another school of thought suggests that the formation of 
fog next to a cold surface is not a direct consequence of the cooling, but rather an effect 
of turbulent mixing of nearly saturated eddies at different temperatures. If a warm air 
parcel advectively flows over a cold surface and does not drive away the cold air 
originally on underlying the surface, but horizontally mixes with the cold air, the 
advection cooling is accomplished by the mixing of warm and cold air. When cold and 
warm air parcels advectively flow into each other, the warm air often vertically mixes 
with the cold air below, and heat (also water vapor) is transported downward due to 
turbulent exchange. Still another hypothesis, based on modeling and measurements of 
fog, states that a virtual cessation of turbulent mixing is necessary before fog can form 
and persist by the radiational cooling of droplets, i.e., turbulent mixing may inhibit fog 
formation at the surface, especially in light winds (Gerber 1981; Turton and Brown 
1987). While other observations have suggested that cessation of turbulent mixing may 
not be required; it is generally accepted that once fog forms, radiational cooling of the 
droplets causes a development and thickening of the fog and that turbulence and 
radiational cooling interact in controlling fog. For example, Turton and Brown (1987) 
suggest that the development of mature fog is dominated by radiative cooling from the 
fog top. Clearly, however, there is still considerable uncertainty as to the processes 
responsible for fog formation (Welch, Ravichandran et al. 1986). 

Advection fog is very dense fog that typically only dissipates when the wind 
direction shifts, or when the fog is advected over warmer land. Over the ocean, 
dissipation occurs when the wind advects the fog across the horizontal SST gradients. 
When warmer water is reached, the increased buoyancy flux often pushes air parcels up 
through the inversion. This causes the entrainment of drier air from aloft to break up the 
fog layer, causing the fog to lift, form stratus, and eventually disappear. 
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Light wind speeds are considered necessary for fog formation and, in fact, fog is 
most often associated with lower wind speed events over land. However, higher wind 
speeds can coexist within advection fog because the relatively smooth sea surface 
generates less friction, resulting in less shear generated turbulence. This type of fog often 
occurs at latitude of about 42° N in the vicinity of the northern wall of the Gulf Stream. 
Once the fog forms over the continental shelf, radiation is much less important in 
controlling this fog than SST gradients. In coastal waters, sea currents act on the 
transportation of waters with different temperatures, thus determining the distribution and 
variation of the temperatures and their gradients. So, they provide the necessary 
conditions for the formation, continuation and dissipation of advection fogs. Therefore 
the primary hydrological element directly related to advection fog should be the sea water 
temperature, especially the SST and its gradient (Binhua 1985; Lundquist 2000). 

The challenges presented by fog to modelers are clearly evident. In summary, the 
strong atmospheric stability associated with fog typically has lighter winds, shorter 
periods of oscillation, and greater likelihood of patchy fog. Dense fog is composed of 
highly structured and sharply defined local regions (Welch and Wielicki 1986). The 
constant flux layer may only be a few meters deep in fogs, with patchy regions of 
turbulence above this layer. Due to the scarcity of turbulence measurements in fogs, 
there is great uncertainty as to the proper parameterization of various coefficients 
required in the higher-order closure formulations (Welch, Ravichandran et al. 1986). In 
addition, boundary layer models usually place the first grid point in air several meters 
above the surface. However, it is believed that important processes in fog formation take 
place below one meter and it is necessary to model them explicitly. Additionally, most 
models are currently unable to accurately reproduce the light winds observed near the 
surface which are essential for fog formation and continuation. This defect is a common 
characteristic of boundary layer models when applied to stable conditions and is believed 
to be due to the finite size of the roughness elements (Turton and Brown 1987). 
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1.2.5 COAMPS 


Models such as COAMPS that rely on flux parameterizations based on MABL 
experiments conducted primarily in unstable conditions (Fairall, Bradley et al. 2003) 
have difficulties predicting fog in stable conditions in the marine environment. As a 
result, the parameterizations for stable boundary layers, particularly when fog is present, 
are very questionable. The inadequate vertical and horizontal resolution is another issue 
which makes model forecasting difficult. For fluxes, COAMPS uses the MOS based 
parameterization from Louis (1979), which is a traditional one and a half-order 
turbulence closure model. This model predicts turbulent kinetic energy and determines a 
mixing length to derive a turbulent exchange coefficient which will be discussed in 
chapter 3. According to one of the COAMPS modelers (Wang 2004), this 
parameterization gives virtually the same results as TC 3.0 over the ocean in stable 
conditions. At present, COAMPS has no particular scheme to forecast fog, although the 
phase change effect in turbulence is included in the TKE equation. Grid-scale liquid 
water (mass and droplet number) is predicted using a two-moment scheme, and fog 
formation occurs when super saturation exists at the first level (10 meters). This is 
common for boundary layer models to place the first grid point several meters above the 
surface. However, as previously mentioned, the disadvantage is that the important 
processes in fog formation take place below one meter (Wang 2004; Binhua 1985; Turton 
and Brown 1987). 

The coupling of CO AMPS is still a work in progress. Presently, the boundary 
conditions for COAMPS are provided by another model, Navy Operational Global 
Atmospheric Prediction System (NOGAPS), which initializes COAMPS at 0000Z and 
1200Z. SST is provided by satellite derived estimates and then COAMPS makes its own 
SST analysis at the surface every time it runs using optimum interpolation techniques. 
With respect to performance, COAMPS has a tendency to overestimate air temperature, 
potential temperature, latent and specific heat fluxes in amplitude, but captures variability 
reasonably well. It underestimates longwave radiation, which compensates for the 
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overestimation of the heat fluxes in the heat budget. This thesis will attempt to aid in the 
identification of some of the physical processes that complicate forecasting for mesoscale 
prediction models such as COAMPS in its current state and provide further insight 
towards the development of improved parameterizations for these processes and the 
ultimate coupling of CO AMPS to an oceanic prediction model. 
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2. EXPERIMENT 

2.1 Data Collection Site 

The experiment was conducted predominantly in coastal waters south of Martha's 
Vineyard (figure 2.1) with an additional component located on Nantucket Island. 

Figure 2.1 CBLAST Data Collection Site 



Aircraft Measurements 


• H2O Sensor 

• rf)? Sensor 



T1 Connection 


The approximately 25 kilometer long southern shoreline of Martha's Vineyard is nearly 
straight with homogenous alongshore topography to the west of Wasque Shoals. The 
shoreline faces the predominant southwesterly winds and seas from the open ocean. 
Therefore, it is an ideal site for investigating MABL physical processes with minimal 
land based influences during extensive periods of onshore winds (Austin, Edson et al. 
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2002). Instrumentation assets deployed during CBLAST included the MVCO, Nantucket 
Island site, two aircraft, the vessel Nobska, the fixed offshore ASIT, five heavy moorings 
(three with full meteorological packages), nine of ten light moorings (one was lost), and 
seafloor based instrumentation. 

The combination of all of these assets provided a time series from seven spatially 
separated locations with meteorological data and eighteen locations with oceanographic 
data. Figure 2.2 illustrates some of the instrument locations in the CBLAST data 
collection region: 


Figure 2.2 CBLAST Offshore Array 



The point measurements are being combined with the spatial surveys from the two 
aircraft and the Nobska to investigate the processes that exchange momentum, heat, and 
mass across and within the coupled boundary layer (CBL) (Edson, McGillis et al. 2003). 
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2.2 Martha's Vineyard Coastal Observatory 


The MVCO includes a small shore lab located 1.5 kilometers inland and a 10 
meter meteorological mast located 50 meters inshore of the shoreline and just behind the 
present location of the dunes. The mast rises approximately 8 meters above the dunes 
and 13 meters above mean sea level. There is also a subsurface node mounted on the 
seafloor in 12 meters water depth, 1.5 kilometers offshore. The node refers to electrical 
components that power the instruments and telemeter the data to shore. The data from 
the meteorological mast, undersea node and associated instrumentation are connected 
directly to the shore lab via an embedded electro-optic power cable. The shore lab is 
connected to WHOI via a T-l data line. 

The site for the meteorological mast is particularly attractive because the gently 
sloping topography at the beach allows sensor deployment above the internal boundary 
layer and most of the flow distortion induced by the shoreline transition from ocean to 
beach and dunes. This places the fast response instruments in marine air for onshore 
winds. The sensors include a 3 axis ultrasonic anemometer which also provides fast 
response temperature measurements derived from its sound speed measurements and an 
IR hygrometer/ C0 2 sensor. Additional sensors measure the mean wind speed and 
direction, relative humidity, temperature, and pressure. The mast also supports a camera 
to visually monitor the cloud coverage and surface wave conditions in its field of view. 

A 10 meter mast extending above the laboratory holds sensors to measure solar and 
infrared radiation, rainfall rate, temperature, humidity, wind speed, and direction. The 
first phase of the MVCO which included the lab, meteorological mast, and offshore node, 
became operational by June 2001 (Edson and McGillis 2003; Austin, Edson et al. 2002). 


22 



2.3 Air-Sea Interaction Tower 


The ASIT is located 3 kilometers due south of Edgartown Great Pond (see figures 
2.1 and 2.2) and spans the water column at a depth of 15 meters and to a height above the 
sea surface of 22 meters (figure 2.3). The tower is connected directly to shore using fiber 
optic cable. The ASIT was completed late in the summer of 2002 and was outfitted with 
an electronics node and directly connected to MVCO in the fall of 2002 to provide data 
transmission and power directly from shore. The ASIT was instrumented starting in the 
spring of 2003 in preparation of the IOP. 


Figure 2.3 Air-Sea Interaction Tower 
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The 22 meter section of the tower in the MABL was equipped with a fixed mast 
that held fast response sonic anemometer/thermo meters paired with IR hygrometers at the 
lowest three levels. The combination of the sonic anemometer, sonic thermometer and 
infrared hygrometer provide DC estimates of the momentum, TKE, SV, sensible heat, 
and latent heat fluxes. 

In addition, sensors were deployed which are capable of measuring static pressure 
fluctuations at two of those levels to estimate the pressure flux. The dissipation rates of 
TKE and SV were computed from these sensors using inertial sub range estimates. Fixed 
sensors capable of measuring the mean profiles were deployed within the vertical array. 
Additionally, a profiling mast was deployed toward the end of the IOP to measure mean 
profiles to go with the turbulent fluxes through mid-October. A distinct advantage of this 
arrangement is that it provides the ability to measure all the terms of the one dimensional 
TKE and SV budgets: 


—u w 


- dU 

dz 


~r~,dV g 1 dwp dwe 

dz 0 V p a dz dz 


= £ 


(2.3.1) 


“7^30 1 dw'O' 2 

dz 2 dz 


—~■ dQ 1 dwq' 2 
w q -—— = N „ 


(2.3.2) 


(2.3.3) 


dz 2 dz 

u + v' + w~) is the TKE; u , v, and w are the three velocity component 


fluctuations; u'w, v'w, w'0', w'q', w'p', w'e, w'0 2 and w q 2 are the kinematic forms of the 
along wind momentum, crosswind momentum, sensible heat, latent heat, pressure, TKE 
fluxes, temperature variance, and humidity variance respectively; p a is the density of dry 

air; £, N e , and N are the TKE, temperature variance, and humidity variance dissipation 

respectively; and and ^ are the mean vertical gradients of velocity, 

dz dz dz dz 

potential temperature, and humidity. The first two terms on the left side of (2.3.1) 
represent the generation of mechanical turbulence through shear, while the third term 
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represents the production (suppression) of turbulence through convection (stratification). 
The fourth and fifth terms neither produce nor consume TKE, instead they act to 
redistribute TKE within the MABL through pressure and energy transport. The first 
terms on the left side of (2.3.2) and (2.3.3) represent production of SV, the second terms 
act to redistribute the variance (Edson and Fairall 1998). 

The DC fluxes of sensible heat, latent heat, and stress are expressed: 


Qh = Pa C pa W ' e ' 

Qe = Pa L e W ' C l 


Pa 


iuw + jv w 


(2.3.4) 


where c p is the specific heat at constant pressure, L e is the latent heat of vaporization of 
water, Q h is the DC sensible heat flux, Q e is the DC LHF, and f is the DC stress vector, 

and w'O', w q, and u'w are the heat, moisture and momentum fluxes respectively 
(Fairall, Bradley et al. 1996). 

The profile measurements from this setup will also be used to provide an in situ 
calibration of the fixed sensors on the rest of the ASIT. Additional measurements above 
the sea surface included instantaneous wave height estimates from laser and microwave 
altimeters; shortwave radiation; longwave radiation; and upwelling brightness 
temperature from radiometers; three more levels with sonic anemometers; a rain gauge; 
and multiple levels with relative humidity and temperature sensors. The data collected 
for this thesis relied primarily on these meteorological instruments located on the ASIT. 
More information on these instruments is listed in table 2.1 where T air is air temperature, 


T v is virtual temperature, RH is relative humidity, U is horizontal wind speed, P is 


pressure, and q air is specific humidity (Edson, McGillis et al. 2003; Edson and McGillis 
2003). 
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Table 2.1 Air-Sea Interaction Tower Meteorological Instruments 


Instrument 

Parameters 

Approximate Height 

(m) 

Sampling Rate (Hz) 

Vaisala HMP45 

T air > q air 

4, 6, 13, 16 

1 

Vaisala PTU 

T ir ,RH, q air ,P 

8, 18 

1 

CSAT3 

U, T v , TKE, Fluxes 

4, 6, 8, 16 

20 

Licor LI7000 

fiair > FluXeS 

13 

20 

Licor LI7500 

fiair > FluXeS 

4, 6,8 

20 

Gill R2A 

U, T v , TKE, Fluxes 

13 

20 

Gill R3A 

U, T v , TKE, Fluxes 

18 

20 

Met3A 

P, Fluxes 

6, 8 

8 

Eppley PSP and PIR 

(Upward) 

Downwelling Solar 

and IR Radiative 

Fluxes 

22 

1 

Wintronics 

Radiometer 

(Downward) 

SST, Upwelling IR 

Radiative Fluxes 

11 

1 

Rain gauge 

Rain 

11 

1 

Riegl Laser 

altimeter 

Wave Height 

Spectra 

11 

20 

Microwave sensor 

Wave Height 

Spectra 

10 

20 

Profiling package 

U, T , RH, q 

9 air 9 9 ~ air 

3-14 

1 
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2.4 Synoptic Situation 


To assist in the interpretation of the measurements, the synoptic meteorological 
situation is investigated during the period of interest from August 20-27, 2003. The 
synoptic conditions are shown in figures 2.4a and b represented by various time series for 
this period measured at the ASIT. These time series are the sensible and latent heat 
fluxes were taken from the CSAT3 sonic anemometers and LI-7500 in fr ared hygrometer 
at 6 meters; temperature, specific and relative humidity were taken from the Vaisala 
HMP 45 probe at 6 meters; specific humidity, temperature (SST) of the sea surface, and 
upwelling IR were estimated from the downward facing Wintronics radiometer deployed 
at 11 meters; downwelling solar and infrared radiation measured from the upward facing 
Eppley PSP and PIR at 22 meters; wind direction (WD1-6) and horizontal wind speed 
(Ul-6) where 1-6 annotate the heights 4, 6, 8, 13, 16, and 18 meters (see table 2.1 and 
figure 2.3 for instrument locations and sampling rates). Significant atmospheric events 
are annotated for correlation with surface meteorological plots and periods are numbered 
for a piecewise discussion. Note that the yeardays are marked in Greenwich Mean Time 
(GMT). The gaps in the wind data of figure 2.4b are a result of dropouts in the sonic data 
during periods of heaviest fog and rain. Periods when fog was present are annotated with 
black stars at the bottom of the composite plots. The criteria for fog at the ASIT that was 
used is as follows: 

\ T air- T bb \^ 6 

q>q s (2.4.1) 

RH > 95% 


where q and q s are the specific humidities for air at 6 meters and sea surface 
respectively; RH is the relative humidity measured at 6 meters; T ajr is measured at 6 
meters and T bb is the blackbody temperature defined as: 


T = 

± bb 


IR 


down 


V ”sb ) 


(2.4.2) 
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Figure 2.4a: Time series showing synoptic variability in sensible heat flux, latent heat 
flux (top panel); air-sea temperature and humidity differences (middle panel); and 
relative humidity (bottom panel). The lines denote significant synoptic events and the 
asterisks denote foggy conditions. 
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Figure 2.4b: As in figure 2.4a but for downwelling IR and solar radiation (top panel); 
wind direction (middle panel); and wind speed (bottom panel). 
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Figure 2.5 Parameters used to identify fog 



where <J sb =5.67x10 8 is the Stefan-Boltzmann constant and IR down is the downwelling 
IR radiative flux measured at 22 meters. This criterion was achieved by matching kn own 
occurrences of fog based on visual observations, to the data. The standard measurement 
level of 6 meters was chosen to minimize distortion by the platform and the expectation 
that it was above the WBL for scalars. A plot of the criteria for fog formation is shown in 
figure 2.5. Temporal means of the measured parameters for each of the periods at a 
height of 6 meters are listed in table 2.1 at the end of this chapter. The data collection 
area consisting of the ASIT and MVCO will hereby be referred to as the data region. 

Period 1 (figures 2.6a-d) : This period was dominated by a high pressure system 
over the data region in mostly clear skies with periods of haze and light fog. It was warm 
and humid with light winds from the west-southwest. There was initially a weak trough 
off the coast of New England on the morning of the 20 th which moved offshore as the day 
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Figure 2.6a 20 Aug 0800L 


Figure 2.6b 20 Aug 2000L 




Figure 2.6c 21 Aug 0800L 


Figure 2.6d 21 Aug 2000L 




progressed. No local precipitation was produced by this trough. A stable boundary layer 
exists in this period and persists throughout the periods of interest. 

Period 2 (figures 2.7a,b) : During this period, the data region was in the warm 
sector of an approaching low pressure system resulting in warm temperatures with mostly 
hazy and foggy conditions throughout the period. Winds were light and from the 
southwest. The specific humidity of the sea was lower than that of the air and both the 
latent and specific heat fluxes were negative during this period. 
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Period 3 (figures 2.8a,b) : The passage of a cold front over the data region was 
evident by a pre-frontal shift of winds to the northwest followed by a post-frontal shift of 
winds to the southwest and the advection of drier air into the region. Skies went from 
partly cloudy to clear as the front progressed eastward with foggy conditions in the early 
morning hours of the 23 ld . Winds were light and predominantly from the west-northwest. 
There was a slightly positive LHF during this period. 


Figure 2.8a 23 Aug 0800L Figure 2.8b 23 Aug 2000L 



Period 4 (figures 2.9a-c) : A high pressure system of relatively cool and very dry 
continental air from Canada resulted in clear, bright skies and a positive LHF. The 
specific humidity of the air was as much as 8.5 g/kg lower than that of the sea. The 
average specific heat flux (SHF) given in table 2.2 indicates stable conditions. However, 


32 










as shown by the time series in figure 2.4a, the period briefly displayed an interval of 
positive SHF and unstable conditions. Wind speed was slightly higher at around 10 m/s 
on the afternoon of the 24 th but tapered off as the wind direction shifted from the north to 
the southwest as a new frontal system approached from the west. 


Figure 2.9a 24 Aug 0800L Figure 2.9b 24 Aug 2000L 



Figure 2.9c 25 Aug 0800L 



Period 5 (figures 2.10a,b) : This period began with data region in the warm sector 
of another approaching low pressure system. Skies went from clear to mostly cloudy by 
the end of the period with the similar characteristic northerly and subsequent southerly 
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shift in the winds as the front passed. This front was not as strong as the previous front 
so the changes in temperature and specific humidity were not as pronounced. 
Nonetheless, conditions were sufficient to meet the criteria for fog formation in the early 
morning of the 26 th . 


Figure 2.10a 25 Aug 2000L Figure 2.10b 26 Aug 0800L 



Period 6 (figure 2.11) : During this period the data region was between the cold 
front which had recently passed over the region and the warm sector of another low 
pressure system to the north. Skies were partly cloudy with light winds from the south- 
southwest. Temperatures were warm with moderate humidity. There was a slightly 
positive LFIF during this period. 

Figure 2.11 26 Aug 2000L 
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Period 7 (figure 2.12) : During this period a trough passed over the data region. 
Winds shifted from the southwest to the west-northwest and were still relatively light. 
Conditions for fog were met late morning on the 27 th and the rest of the period was 
characterized by partly cloudy skies. 

Figure 2.12 27 Aug 0800L 



Period 8 (figure 2.13) : This brief period at the end of the case study is 
characterized by an approaching cold front. Winds were between 5-10 m/s and from the 
west. Temperatures were warm with relatively mild humidity. 


Figure 2.13 27 Aug 2000L 



35 










Table 2.2 Period Means 


Parameter 

1 

2 

3 

4 

5 

6 

7 

8 

T air (C) 

21.86 

22.63 

22.47 

19.79 

20.65 

20.79 

20.82 

21.91 

q ai r(g /k g) 

14.76 

16.41 

13.74 

7.84 

13.10 

13.20 

14.86 

13.79 

SST (C) 

20.77 

20.66 

20.52 

19.13 

19.48 

20.17 

19.61 

19.62 

q sea (g /k g) 

14.89 

14.89 

14.81 

13.50 

13.81 

14.45 

13.99 

13.98 

LHF(W/m 2 ) 

4.06 

-12.31 

16.57 

106.73 

9.45 

16.24 

-4.12 

3.87 

SHF (W/m 2 ) 

-6.60 

-11.91 

-8.35 

-1.04 

-5.69 

-1.44 

-7.90 

-13.61 

RH (%) 

90.77 

95.71 

81.01 

54.73 

86.90 

86.28 

96.51 

84.16 

U (m/s) 

4.84 

5.24 

4.36 

5.37 

4.37 

4.34 

4.89 

5.94 

WD (°) 

249.84 

225.04 

283.23 

262.18 

231.88 

194.17 

217.66 

254.35 
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3. THEORY 


3.1 Monin-Obukhov Similarity 

Monin-Obukhov Similarity (MOS) is a statistical tool that was developed to study 
atmospheric turbulence in the surface layer over land and is sometimes referred to as 
surface-layer similarity (Monin and Obukhov 1954). The surface layer is defined as that 
part of the boundary layer where the fluxes vary by less than 10% and therefore MOS 
assumes a nearly constant flux layer. MOS also assumes that the statistics are stationary 
and horizontally homogenous. Within the boundary layer, and near the surface, turbulent 
energy is produced mainly by mechanical working of the stresses on the mean velocity 
gradient and higher up, principally by buoyant motions (Bradley, Coppin et al. 1991). 

The fluxes associated with these two production mechanisms are represented by the 
following: 

— 2 

— =k(° )f (3.1-1) 

Pa J 

Vei = -^- + 0.6IT, -^ = F.(0) + 0.6irF (0) (3.1.2) 

5 Pa C p Pa L e 

where f s is the stress vector at the surface (i.e., the surface stress), w, is the friction 
velocity, Q h is the surface value of the sensible heat flux, Q e is the surface value of the 
LHF, T s is the sea surface temperature (SST) and F u ,F e , and F q are the kinematic forms 

of the momentum, sensible heat, and latent heat fluxes, at the surface respectively. 

The friction velocity and heat fluxes are combined with the buoyancy parameter, 

a 

-2-, where g is the acceleration due to gravity and 0 is the mean virtual potential 
0 

temperature; and the height above the surface, z , to form a set of scaling parameters for 
velocity, temperature, humidity and an additional length scale. The scaling parameters 
can be defined as: 
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X, 


(3.1.3) 


.. MU 


u t 


where x* is the respective scaling parameter forx = q,6,u . The additional length scale 
that results from a Buckingham-Pi analysis (Sabersky, Hauptmann et al. 1998) is known 
as the Obukhov length, L . The Obukhov length is proportional to the height at which 
buoyant production equals mechanical (shear) production of turbulence in a constant flux 
layer (2.3.1), and is expressed as: 


© v w* _ 0 v w* 
g tew'd', g>c(0 t + 0.61 T s q ,) 


(3.1.4) 


The Obukhov length and height z are then combined to form the dimensionless 
parameter, £ , called the surface layer scaling parameter: 


z _ ~Kzg(w'd[) 

L Oul 


(3.1.5) 


The sign of the surface layer scaling parameter relates to static stability: negative implies 
unstable conditions, and positive implies stable conditions. 

The basis for MOS is that flows with similar ratios of convective to mechanical 
generation of turbulence at a given height (i.e similar £ or Richardson numbers) should 
have similar statistical properties after normalization by the appropriate scaling 
parameters. Specifically, the similarity hypothesis states that normalization of surface 
layer variables by the appropriate scaling parameter(s) should be universal functions of 
the surface layer scaling parameter £ . For example, using these scaling parameters, the 


TKE budget (2.3.1) can be normalized by — to obtain: 

u t 


Ut 

where the (j) functions represent the dimensionless form of the terms in (2.3.1). In a 
constant flux layer with no turning of the wind, the dimensionless gradient for velocity 
<j) u becomes: 
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(3.1.7) 


t{C) 


Kz dU 
lit dz 


Likewise, the dimensionless gradients of specific humidity and potential temperature are; 


MO 

MC) 


Kz dQ 
cjt dz 
Kz 30 
0,. dz 


(3.1.8) 


where U,Q, and 0 are the mean values. As such MOS predicts that these non- 
dimensional gradients of velocity, humidity, and temperature are universal functions of 
atmospheric stability, £. 

These relationships are commonly used to relate the fluxes to their respective 
profiles. For example, the combination of (3.1.3) with (3.1.7) and (3.1.8) results in the 
flux parameterization: 


F x = w'x = -u. t Xt 


U tKz dX _ dX 
</>AC) dz 1 dz 


(3.1.9) 


where again, x = q,6,u , and K x is the respective turbulent exchange coefficient. This is 
a first order, small eddy closure technique implying a down-gradient transport. These 
scaling laws are expected to hold and the derived parameterizations are expected to be 
universal as long as the assumptions that govern MOS laws are valid, i.e., a combination 
of mechanical and thermal forcing drive the turbulent exchange, the scaling parameters 
are independent of height in the surface layer, and the turbulence statistics are stationary 
and horizontally homogenous (Edson, Zappa et al. 2002; Stull 1988). 

In the MABL, one can think of several situations where the assumptions 
governing MOS could become invalid. For example, very near the ocean surface, we 
expect to encounter a wave boundary layer (WBL) as described in section 1.2.2, where 
U(t)=U +u'(t)+u(t) and w (/) represents the wave induced fluctuations in the velocity 

field. In such a flow field it is generally assumed that MOS is not valid since an 
additional forcing mechanism influences the near surface flow, i.e., the wave induced 
flow (Sullivan 2004; Edson, Zappa et al. 2002). 
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The primary question to be addressed in this thesis is whether or not MOS is valid 
in stable conditions. In the stable regime, the stratification begins to restrict the 
production of TKE through shear by limiting the size of the energy containing eddies (i.e. 
their velocity fluctuations). This is a result of the restoring forces that limit the 
displacement of the parcels from their equilibrium position. In extremely stable 
conditions, the size of the eddies are completely limited by the stability and they become 
unaware of their distance from the surface. The scaling becomes height independent 
under these conditions and the Monin-Obukhov (MO) length becomes the only length 
scale. As a result, we often refer to such conditions as z-less stratification (Welch, 
Ravichandran et al. 1986; Wyngaard, Busch et al. 1973). The weak and intermittent 
turbulence often observed in stable conditions, particularly in low wind conditions, leads 
to non-stationarity and inhomogeneity, which strongly impacts our ability to model the 
turbulent mixing of heat and water vapor. The processes controlling fog within the 
surface layer further complicate the parameterization of these fluxes (Gerber 1981). 

3.2 TOGA COARE 3.0 Parameterization 

The TOGA COARE 3.0 implementation of the bulk aerodynamic formula 
provides a starting point for this investigation. As stated earlier, TC 3.0 has been 
validated in several field programs characterized by unstable to near neutral conditions. 

As opposed to the DC fluxes described in section 2.3, the bulk aerodynamic 
formula parameterize the kinematic fluxes in terms of the more easily measured mean or 
bulk quantities and are expressed: 

Q,=p.c ra c„s(T,-e) 

Q,=P.L,C E S(q,-q) (3.2.1) 

r, = p„C D S(u s ,-u,) 

where 0 , q,u t and S are the average potential temperature, specific humidity, horizontal 
wind velocity in the i th direction, and instantaneous wind speed, respectively, at some 
height z ; ,v denotes their surface values and C H ,C E , and C D are the transfer coefficients 
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for sensible heat (i.e. the Stanton number), latent heat (i.e. the Dalton number), and 
momentum (i.e. the drag coefficient) respectively. The interfacial value of the water 
vapor mixing ratio, q s , is: 

(3.2.2) 

where the value of 0.98 multiplying the saturation specific humidity of the SST [ q sat ( T s ) ] 

takes into account the reduction in vapor pressure caused by a typical salinity of 34 parts 
per thousand. In this analysis, the potential temperature and specific humidity at 
reference height z r are determined from: 


6 = T + 0.0098z r 

9=m,jT) 


(3.2.3) 


where RH is the relative humidity and T is the air temperature and 0.0098 K/m is the 
adiabatic lapse rate. 

The use of S in equations (3.2.1) is an important point explained in the following 
excerpt from Edson (2003). As buoyant production begins to dominate the generation of 
turbulence in a very unstable (convective) atmospheric boundary layer, the mean wind 

vector, U , and the mean shear, , approach zero. However, the variance of the 

oz 

velocity components and instantaneous shear remain finite due to the convective motion 
that drives eddies capable of transporting heat and momentum. The size of these eddies 
scales with the boundary layer depth and the strength of the buoyancy flux, and define a 
convective velocity scale given by: 


f 


w» 


— w'6' v z 

\ d v 




(3.2.4) 


J 


The TC 3.0 bulk algorithm attempts to account for the heat and momentum exchange 
driven by these convective eddies by modifying the traditional definition of the bulk 
formula to include the mean instantaneous wind speed: 





(3.2.5) 
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where the subscript r denotes that the values are relative to the sea surface. 

Unfortunately, it is not always possible to accurately measure the mean wind 
speed over the ocean due to platform motion. For example, the mean wind speed 
measured from a buoy mounted cup anemometer will include fluctuations from both the 
wind and wave induced motions. For this reason, the mean wind components (i.e., vector 
averaged winds) are normally measured on ships and buoys because the sinusoidal nature 
of the wave induced motions tend to average out. The following definition is then used 
to estimate the wind speed from the mean wind components: 


s = { u l + u l+ u l 


\ 2 1—2 , —2 . 2 \2 

j =[u rr +u„. +w„ 


ry 


(3.2.6) 


where w g is known as the gustiness and provides an estimate of the wind speed 

measurements of the mean wind vector. The gustiness is related to the convective 
velocity by: 

w g =j3w t (3.2.7) 

where j3 is a numerical constant known as the gustiness parameter with a value of order 
1. Fortunately for this study, the fixed sensors on the ASIT allow direct measurement of 
S. 


The parameterization of the transfer coefficients begins with commonly used 
forms of the diabatic profiles found through integration of (3.1.7). This provides the 
mean quantities at height z : 


U(z) = ^ 
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where the stability parameter y/ x is defined as: 

£ (\-</> (0) 

y,(0 = J r dC 

0 £ 
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and z o e are the roughness lengths for momentum, heat, and moisture respectively. The 

roughness lengths are defined as the height where the extrapolation of the log-z portion 
of the respective profile intersects the surface value. The most commonly used forms of 
the dimensionless gradient functions are known as the Businger-Dyer formulations and 
are given by: 


&=i+or ] 

^=^=i+c 2 <rj 
A=(i+Q0 4 
^ = ^=(i+Q0 4 


0 < ^ < 0.2 


-1 < ^ < 0 


(3.2.10) 


where C # are the numerical coefficients found from field experiments. Decades of 
overland experiments have resulted in a range of values. However, commonly used 
values are given by C l = C 2 = 5 and C 3 = C 4 = 16, which can be combined with (3.2.9) to 
define the stability parameters: 


V u {0 = Ve{0 = V q {?) = 5?} stable 


(3.2.11) 


MO = 21n 


1 +( 1 - 160 " 


+ ln 


1 + ( 1 - 160 " 


n 


-2tan" 1 (l-16^'p+- 


MO = MO = 21n 


1 + ( 1 - 160 " 


unstable (3.2.12) 


(Paulson 1970; Large and Pond 1982): 

The TC 3.0 algorithm partitions the transfer coefficients into individual profile 
components: 

1 1 

c — r 2 r 2 

C E = c\cl (3.2.13) 

1 1 

c — c 2 c 2 

The combination of (3.2.13) with (3.2.1) and (3.2.8) yields the following semi-empirical 
forms of these components: 
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(3.2.14) 


where a accounts for the difference in scalar and velocity von Karman constants. The 
neutral transfer coefficients defined by £ = 0 and lf/ x = 0 are related to these transfer 


coefficients through (Large and Pond 1982): 
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(3.2.15) 


The MOS scalar scaling parameters can be computed independently from the velocity 
measurements: 


6t = ~Cg (T s - 9) 

V* = ~4 (? s - 9) 


(3.2.16) 


The ability to compute the scalar components separately is advantageous because the 
drag coefficient is sensitive to both sea state and wave age, while the scalar coefficients 


44 



may be influenced by additional processes such as wave breaking and heat exchange 
from evaporating sea spray (Edson, Zappa et al. 2002). 

Because the transfer coefficients are functions of height and stability it is common 
practice to adopt a reference height of 10m and adjust profiles to their form at neutral 
stability to compare results from different elevations. The expression "neutral transfer 
coefficient" is a classic contradiction in terms, since neutral stability implies zero heat 
flux at the surface and a nonexistent potential temp gradient (Bradley, Coppin et al. 

1991). In addition to neutral transfer coefficients, the neutral wind speed is found from: 


U„ = U}S h + ¥ "^ u * = u) (U; + w 2 g f + 


(3.2.17) 


K ' K 

The velocity roughness length z o is often related to the viscous sub layer at low 


winds and the physical roughness of the surface at higher winds. The scalar roughness 
lengths are related to the thermal sub layer. It has proven convenient to characterize the 
surface and the flow regime by the roughness Reynolds number: 


R = 


u*z 


V 


(3.2.18) 


As the wind speed decreases, laboratory experiments have shown that ^.approaches a 
constant value of about 0.11 and the relationship between roughness and stress is given 
by: 


0.1 In 


(3.2.19) 


The validity of this smooth flow relationship at low winds over a malleable ocean surface 
(i.e, its surface characteristics are also governed by surface tension) remains a matter of 
some debate for field applications. 

The relationship between oceanic roughness and stress for rough flow is: 


z„ = - 


ecu , 
g 


(3.2.20) 


where eris the Charnock constant which has values between 0.010 and 0.035 and is 
linked to gross characterizations of the sea state such as wave age or slope of the 
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dominant wavelength (from the peak of the gravity wave spectrum). Combining (3.2.19) 
and (3.2.20) we get a more comprehensive expression forz o (Smith 1988; Fairall, 
Bradley et al. 1996): 


aul 0.1 \v 

Z o = -+- 

g Ut 


(3.2.21) 


The dependence of the Charnock "variable" as a function of sea state is the focus of many 
ongoing investigations 


3.3 Local Similarity 


In section 3.1, a condition referred to as z-less stratification is described for the 
stable regime as the condition when the buoyancy force begins to restrict the production 
of TKE through shear by limiting the size of the energy containing eddies, i.e. their 
velocity fluctuations. In other words, local similarity (LS) recognizes that turbulence in 
the mid and upper stable MABL may not be in equilibrium with the surface fluxes. This 
means that the local fluxes, shears and stability are more important than the surface 
values. For comparison, the orders of magnitude for MOS scaling factors in unstable 
conditions are: 


L 

z 


w* 


6 1 


q* 


-O ( 1 in - 200/m) 

~0(0.0 1/m - 0.1z,/m) 


TVL 111 

O 0.05— -0.3 — 
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while for LS, the relevant scales are: 

L (z) 

w*( z ) 

0\ (z) 

q*( z ) 


O(0m - 50 m) 

- O(0m/s - 0.3 mls ) 

- O(-2.0°C - 0°C) 

“ 0(-2g / kg - 5 g/kg) 


Dimensionless groups formed from these scales are not functions of height z , although 
the individual scaling variables listed above may vary significantly with height (Stull 
1988). These LS scaling factors are computed as: 

u L = u, (z) 

L _ ® v u L {zf 

g/cw'e'(z) 


W0' v (z) 


(3.3.1) 


wq'{z) 

u L {*) 


where L l is the only length scale. 


In conditions of z-less stratification, the velocity, temperature, and humidity 

X 

gradients are proportional to —, i.e: 

L, 


and therefore: 


dz L, 


KZ dx KZ X, z 

(p =-oc-oc- 

x, dz x, L l L l 


Most experiments show that in stable conditions the following is true: 

A=1+/? i 


(3.3.2) 


(3.3.3) 


(3.3.4) 


which is the same as equations (3.2.10) with j3 = C # (Wyngaard, Busch et al. 1973). 
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All of these relationships assume that the stable MABL is continuously turbulent 
in time and space with no gaps or patches of nonturbulent air. Since real stable MABLs 
have sporadic, patchy turbulence, it is important to recognize that these expressions have 
their limitations (Stull 1988). 

3.4 Enthalpy and Specific Enthalpy 

To determine our best estimates of the DC fluxes for comparison with the TC 3.0 
parameterizations of the scalar transfer coefficients, we examine the enthalpy exchange 
of the system. The following discussion of enthalpy is from Edson (2004) based on the 
work of Businger (1982) and Webb et al. (1980). Enthalpy is a combination of the latent 
and specific heat fluxes and therefore the exchange of enthalpy is of paramount 
significance in representing an estimate of the total energy entering a system. As a result, 
care should be taken to use consistent forms of the sensible and latent heat fluxes to 
ensure that enthalpy is conserved. Through a detailed investigation of the enthalpy 
exchange in stable conditions with and without fog, we hope to gain valuable insights 
that will point to possible sources for error in the TC 3.0 parameterizations. This 
discussion presents the traditional equations for enthalpy and then discusses correction 
factors that may be included in the equations to account for internal energy changes due 
to the presence of water vapor and an additional correction for the presence of liquid 
water droplets. Ideally, the addition of the correction factors and the correction for the 
presence of liquid water droplets should improve the accounting for the complex energy 
exchanges associated with fog in the MABL. 

3.4.1 Traditional Equations 

We begin with the traditional equations defining enthalpy from Emanuel (1995): 

k d =c pd T (3.4.1) 

K=c p J (3.4.2) 
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(3.4.3) 


k =c T = k +L =c T + L 

v pv we pw e 

where k d ,k v ,k w and c d ,c ,c are the specific enthalpy and specific heat at constant 

pressure of dry air, water vapor, and liquid water, respectively; and T is the temperature 
of these constituents. 

Even if no fog or sea spray is present in the MABL, there is still moisture present. 
So we begin with the equation for specific enthalpy of moist air just above the surface as 
defined by Emanuel (1995) using equations (3.4.1), (3.4.2) and (3.4.3): 

k m={ l -‘l) c P d T + c t c p V r ( 3 - 4 - 4 ) 

or 

K=0-9KT+q( Cl J+L.) (3.4.5) 

Using (3.4.5) and the expansion described by Businger (1982) we can derive the surface 
enthalpy flux as: 

pwk nh » pw'T'[( 1 ~q)c pd + qc pv ] + pw q \c p T + L e ] (3.4.6) 

3.4.2 System Enthalpy 

To investigate the enthalpy carried by a parcel of air and how it impacts the 
enthalpy of the system away from the surface, Businger (1982) included a set of arbitrary 
constants in the equations for enthalpy as follows: 


k d= c pj+ b d 

(3.4.7) 

K= c P v T + b v =k w +L e 

(3.4.8) 

K = C P J+K 

(3.4.9) 


where b d ,b v ,b w are constants which are dependant upon the internal energies of these 
constituents. The latent heat of vaporization, L e , is not considered a constant here, but a 
function of temperature, L e (T). 


49 



Using equations (3.4.7), (3.4.8), and (3.4.9) we obtain the following expressions 
for specific enthalpy of moist air: 


k m=( l - ( l) C P d T+( l C pv r + b mX 

(3.4.10) 

or 


K =(!-?) c P d T +<i c pJ + qL e + b m2 

(3.4.11) 

where 


Kx = { l ~q)b d + qb v 

(3.4.12) 

h ,n2={ [ - C l) b J + C lK 

(3.4.13) 

Now an enthalpy flux can be derived using (3.4.11) and (3.4.13) as follows: 

pwk m = pw{\- q)c pd T + pwqc p J + pwqL e + pwb m2 

(3.4.14) 


on which an expansion is performed based on Businger (1982) but modified by Edson 
(2004) to be consistent with the results of Webb et al. (1980) to obtain (see Appendix A): 

pwk m =pw'T'[(\-q)c pd + qc pv ] + wq \c p f + L e +b w ~\ (3.4.15) 


As Businger (1982) describes, the combination of constants with the specific 
humidity in the previous expressions are necessary because the source of the water vapor 
in a parcel of air may have different initial conditions. Since the sea surface is the 
primary source of water vapor in the MABL even in the absence of sea-spray or fog, one 
needs to consider the specific enthalpy between the sea surface and the water vapor at the 
surface. Businger (1982) asserts that this is equivalent to the expression: 

b =-c T (3.4.16) 

where T o is the average sea surface temperature. This closes the expression for the moist 
enthalpy flux in the surface layer, and provides an expression for the enthalpy flux away 
from the surface: 

pwk m ^pw'T'[(l-q)c pd +qc pv ] + j^-^wq'[c pw (T-f o ) + L eo ~\ (3.4.17) 


50 



Alternatively, we can use the relationship L e —L e = [c pr - c ;i i ) (T - T o ) to obtain: 

pwk m ^pw'T\{\-q)c pd +qc p ^ + -^-^wq[c pv (T-T o ) + L e ~\ (3.4.18) 

which is Eq.(17) in Businger (1982) and is consistent with Eq.(39) in Webb et al. (1980). 
As explained in Fairall et al. (1996), the surface temperature in equation (3.4.17) 
represents the heat required to cool the water vapor from its initial temperature T o to the 

air temperature T for unstable conditions. Webb et al. (1980) describe this term by 
noting that the heat imparted to and carried by the air parcel is represented by the change 

in temperature (T -T o j , not the temperature itself. 
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(3.4.24) 


h = - 


pwk m 


U* 


where k m is the specific enthalpy of the sea surface. Additionally, we can define the 
dimensionless enthalpy flux as follows: 


Kz dk m 

__ m_ 

k, dz 


(3.4.25) 


and 


dk„ 




(c pi -qc rl + qc„)^+[(T-T„) Cpw -c^T+L,]^ (3.4.26) 


3z ' ' ,jz 

which was derived from (3.4.11) using b m2 = b w = —c T 0 such that: 


dz 


( x ~q) c P d T + <l\c pw (T ~T 0 ) + L e 


(3.4.27) 


3.4.3 Inclusion of Liquid Water 

In fog, there are liquid water droplets present in the air in addition to water vapor. 
The discussion of enthalpy thus far only accounts for the effects of water vapor on the net 
enthalpy flux and does not consider the effects of liquid water droplets. Therefore, we 
now consider the effects of liquid water on the net enthalpy flux. We start with an 
expression from Frank and Emmitt (1981) for the specific enthalpy of moist air that 
includes liquid water: 

K = ( :1 - q ~ q w ) K + qk v + q w k w (3.4.28) 

or 

K = 0-7- q w )k d +(q + q w ) k w + qL (3.4.29) 

where q w is the specific enthalpy of liquid water. The moist enthalpy is then given by: 

k m = (1 - q ~ q w )c pd T + {q + q w ) cj + qL e + b mi (3.4.30) 

where 

b m 3 = 0 - q ~ 7, )b d +{q + q w )b w (3.4.31) 
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From (3.4.30) the moist enthalpy flux including liquid water is: 


pwk m = pw{\ — q — q w ) c pd T + pw(q + q w )c p J + pwqL e + pwb m3 (3.4.32) 

Using equation (3.4.32) and the expansion given in Appendix B, we obtain an estimate of 
the moist enthalpy flux including liquid water; 


pwk m = pw'T'[(l-q- q w ) c pd + qc pv + q w c pw ] 


+ - 


P 


( l ~q-q w ) 

The sensible and latent heat fluxes then become: 

P 


Qh =pw'T'[(\-q-q w )c pd +qc pv +q w c pw ] + 

P 


Qe 




( 1 -? 

L Vo ( w'q' + w'q ' H ,) 


(wq' + w'q w ) \c pw (T -T 0 ) + L % ](3.4.33) 

c„(f-2;)(3.4.34) 
(3.4.35) 


w'q'+w'q' H , 


g 

where the liquid water content in typical marine fog is q = 0.05 - 0.1 (Wallace and 

m 

Hobbs 1977) and w'q' w is the liquid water flux. 

Equations (3.4.34) and (3.4.35) for the SHF and LHF respectively appear to 
represent a full accounting for water vapor and liquid water in the air. If these equations 
are incorporated into the TC 3.0 parameterizations, one would expect that a more 
accurate estimate of the sensible and latent heat fluxes would be achieved not only for 
stable conditions, but stable conditions with fog as well. The next chapter will 
investigate this premise and attempt to close the error gaps between TC 3.0 and DC 
measurements under these conditions. 
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4. RESULTS AND DISCUSSION 


4.1 Transfer Coefficients 

We now return to the question as to whether MOS is valid for parameterizing 
scalar fluxes in stable conditions with fog. Due to the efforts of scientists through the 
various field experiments listed in section 1.1, we know that MOS works well in unstable 
conditions where turbulence is driven by buoyancy and shear (Edson 2004). In contrast, 
the applicability of MOS in strongly stratified boundary layers remains an open question. 
Our poor understanding of the stratified boundary layers has motivated several recently 
completed overland field experiments to investigate this and other questions (Poulos, 
Blumen et al. 2002). In addition, the effects of fog on the net enthalpy exchange between 
the ocean and the atmosphere are not well known due to the sparseness of data collected 
under these conditions. Fortunately, the data collected during the case study period of 
August 20-27, 2003 contains measurements in stable conditions with fog and should 
allow us to investigate this issue to an unprecedented degree of detail. 

To begin the examination of the data acquired during the case study period, we 
return to the MOS based TC 3.0 bulk equations (3.2.1). By inspection, unknowns in the 
bulk parameterizations are the transfer coefficients C E and C H which can be calculated by 
equations (3.2.13) and (3.2.14). Therefore, the proper parameterization of these 
coefficients is essential to the accuracy by which the TC 3.0 bulk equations can predict 
the SV fluxes. 

4.2 Moisture Corrections 

Figures 4.1 through 4.6 are plots of the neutral values of C E ,C H and C K corrected 
to 10 meters versus the neutral value of the mean horizontal wind component corrected to 
10 meters are shown in the top panels. The lower panels show the^,^, and <j) k versus 
stability with no fog present as calculated with equations (3.2.13), (3.2.14), (3.2.15) and 
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(3.1.7), (3.1.9) respectively. The foggy conditions have been removed using the criteria 
developed in section 3. The direct covariance fluxes used to calculate the transfer 
coefficients and the dimensionless fluxes were calculated using the equations for the 
latent and specific heat fluxes (3.4.19) and (3.4.20) as described in section 3.4.2 by 
Businger (1982). 

The plots show a comparison of the direct covariance values versus the TC 3.0 
values. The accuracy of the TC 3.0 parameterizations for the neutral transfer coefficients 
in both unstable and stable conditions with no fog present can be seen in the upper panels 
of figures 4.1, 4.3, and 4.5. The bottom panels are the means and standard deviations for 

z 

stable conditions only. Since the dimensionless profile functions are functions of—, it 

makes sense to describe the ability of TC 3.0 in parameterizing the fluxes by looking at 
the unstable and stable conditions separately. In unstable conditions with no fog present, 
TC 3.0 performs reasonably well in parameterizing all of the transfer coefficients. The 
Dalton number is nearly prefect in unstable conditions as can be seen in figure 4.1. The 
Stanton number is only slightly overestimated in unstable conditions as can be seen in 
figure 4.3. Figure 4.5 shows that TC 3.0 slightly underestimates the enthalpy coefficient 
in unstable conditions. The behavior of the dimensionless flux parameterizations in 
unstable conditions are good as well. There is very good agreement for (f) q and (f) B in 

figures 4.2 and 4.4. In figure 4.6, we see that TC 3.0 slightly overestimates (f) k . Despite, 

these minor differences, these plots validate the conclusions of recent field studies that 
the TC 3.0 parameterizations work remarkably well in unstable conditions. Additionally, 
ongoing investigations to improve the calibration of the LI-7500 indicate that the 
recalibration will reduce these differences even further. 

In stable conditions with no fog present, TC 3.0 does a fair job of parameterizing 
all of the transfer coefficients, but with a slightly larger difference with direct estimates. 
In stable conditions we can see a slight overestimation of both the Dalton number and the 
Stanton number as in the bottom panels of figures 4.1 and 4.3. From the bottom panel of 
figure 4.5, we see that TC 3.0 shows good agreement for C K in stable conditions. 
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Figure 4.1 Neutral Dalton Number 
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Figure 4.2 Dimensionless Moisture Flux 
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Figure 4.3 Neutral Stanton Number 
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Figure 4.4 Dimensionless Heat Flux 
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Figure 4.5 Neutral Enthalpy Coefficient 
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Figure 4.6 Dimensionless Enthalpy Flux 
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Although upon closer examination one could argue that TC 3.0 slightly overestimates C K . 

A possible source of error in the TC 3.0 parameterizations of C E ,C H , and C K 
could be due to the improper parameterization of </> q ,</> g , and </> k in stable conditions. By 

examining the behavior of the TC 3.0 dimensionless flux parameterizations in stable 
conditions we can clearly see that it is not quite as good as the parameterizations in 
unstable conditions. In figure 4.2, TC 3.0 consistently underestimates (j) q but not to a 

great degree. In figure 4.4, TC 3.0 does fairly well with (/) e in near neutral to low 
stability but appears to degrade by increasingly overestimating as stability increases. In 
figure 4.6, TC 3.0 consistently overestimates (j) k in stable conditions. Since enthalpy is a 
combination of LHF and SHF, one would generally expect a direct correlation between 
the errors in the (j) q and </> 0 parameterizations and the error in the (j) k parameterization. 

Based on the results in figures 4.2, 4.4 and 4.6, it appears that q\ has a greater effect on 
<j) k than does </) q in stable conditions with fog. This implies that the sensible heat flux 

dominates the enthalpy exchange when the flow is highly stratified in this data set. 
Nonetheless, these results suggest that the presence of water vapor does indeed affect the 
enthalpy and therefore the net energy exchange between the ocean and the atmosphere. 
Even with the correction factors added to the equations for the enthalpy, the discrepancy 
still persists. It also appears quite possible that the deficiencies in the TC 3.0 
parameterizations for (j ) q , (f ) f , and <j> k might be possible sources for the deficiencies in the 

TC 3.0 parameterizations for C E ,C H , and C K . This relation is illustrated in equations 

(3.1.9), (3.2.9), (3.2.15), (3.2.14) and (3.2.13). However, other possible sources of error 
may indeed be present. 
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4.3 Liquid Water Corrections and Fog 


The instruments used to measure w q were the LI-7500 open path hygrometers 
listed in table 2.1. According to a design engineer, "the measurement principle and the 
measurement electronics of the LI-7500 should focus on the gaseous phase of water 
vapor. Water droplets in the LI-7500 measurement path will scatter "some" light on an 
absorption phenomena but the net effect on the output from the H 2 0 channel when water 
vapor is super saturated should be minimal, but it is not zero. What that small non-zero 
value is, is difficult to quantify." (Anderson 2004). As a result, the LI-7500 in correlation 
with the sonic anemometers gives the water vapor mass flux w p' v with a small amount of 
error due to the presence of the water droplets. Unfortunately, neither liquid water nor 
the liquid water flux were measured during the experiment. Therefore, for the purposes 

of calculations, we further assumed that w q' w « w q , and ignored the liquid water flux 

in equations (3.4.34) and (3.4.35). To estimate the liquid water content, we used the 
results given by Wallace and Hobbs (1977), as shown in figure 4.7. Wallace and Hobbs 

_ a _ 

(1977) report values of /? v between 0.05 and 0.1 -fy- in fog, which is equivalent to q w 

m 

between 0.04 and 0.8 —. The larger value was then incorporated into equations (3.4.34) 

kg 

and (3.4.35) with our assumption of zero liquid water flux for foggy conditions. 

Figures 4.7 through 4.12 are the same plots as figures 4.1 through 4.6 with the 
data from foggy conditions calculated using equations (3.4.34) and (3.4.35) plotted over 
them. With respect to the previous discussion of the results shown in figures 4.1 through 
4.6, the results of the additional consideration of water droplets in the equations show a 
very slight improvement in C H , but a slight degradation in C E and C K when compared 

with the fog-free conditions. Although the inclusion of q w slightly reduces the 
discrepancy between fog and fog-free conditions, we see that in fog TC 3.0 
underestimates (j) q more than it did in stable conditions without fog. Likewise, in figure 
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Figure 4.7 Neutral Dalton Number with Liquid Water 
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Figure 4.8 Dimensionless Moisture Flux with Liquid Water 
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Figure 4.9 Neutral Stanton Number with Liquid Water 
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Figure 4.10 Dimensionless Heat Flux with Liquid Water 




62 







































Figure 4.11 Neutral Enthalpy Coefficient with Liquid Water 
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Figure 4.12 Dimensionless Enthalpy Flux with Liquid Water 
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Figure 4.13 Fog Droplet Concentration Taken from Wallace and Hobbs, 1977 



4.10, TC 3.0 overestimates (f) ft more than it did in stable conditions without fog. In figure 
4.12, these differences appear to cancel each other in the calculation of (j) k as there is very 
little change from stable conditions without fog. These results indicate that the inclusion 
of q w cannot explain the apparent differences. Further investigations are hampered by the 
lack of liquid water fluxes. 
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5. CONCLUSIONS 


5.1 Discussion 

The goal of this thesis was either to prove the incapability of the MOS based TC 
3.0 parameterizations in their present form to account for the proper representation of the 
scalar fluxes, and to identify other possible mechanisms that may be causing observed 
differences between the data and parameterization. The results presented in thesis 
present strong evidence that MOS based TC 3.0 parameterizations for scalar fluxes, the 
moisture flux in particular, are indeed insufficient in stable conditions with fog. 

This case study, which was a subset of the CBLAST experiment, shows that 
systematic differences between the observations and parameterization persist throughout 
all the corrections applied during the course of the analysis. One possible source of error 
for this deficiency is instrument error. As mentioned in section 4.2.3, the LI-7500 
hygrometer is not expected to measure the liquid water flux. However, the effects of 
liquid water on the instrument reading are relatively unknown. To identify and quantify 
this error will require more specialized measurements or techniques. Additionally, it is 
not entirely certain that there was no trace of condensation on the infrared windows of the 
hygrometers when fog was present. Future work will incorporate the internal instrument 
diagnostic logged during the experiment. The use of these diagnostics will help to 
remove data when condensation is a problem. 

Our analysis has shown that the distinct difference between TC 3.0 and the 
observation persists in stable conditions even when the foggy cases are removed. 
Therefore, a more intriguing explanation of this discrepancy is the possibility that there 
are still microscale physical processes present in fog that we still have yet to discover. 
Recall that we have measured all the terms in the 1-D scalar variance budget. This 
investigation will therefore continue by using this budget to determine, e.g., the relative 
importance of the individual terms and how they compare with previous field studies 
(Edson and Fairall 1998). If the most important processes in fog formation and 
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development occur below one meter, perhaps more measurements closer to the sea 
surface will provide further insights. This, of course, becomes quite challenging due to 
swell and wave action. 

5.2 Future Work 

The results of this case study have shown that there is much work yet to be done 
towards understanding the complex physical processes taking place in stable conditions 
with fog. The ultimate goal of future work in this field is the improvement of the TC 3.0 
parameterizations under these conditions. One of the first steps to be taken toward this 
goal would be a closer look at how the instruments were calibrated. As any field scientist 
knows, an essential factor in any field experiment is proper instrument calibration. If the 
instruments used to measure the moisture fluxes during this experiment were not properly 
calibrated, this would certainly introduce errors. More precise calibration techniques will 
be applied to the data collected during this case study and throughout all of the data 
collected during CBLAST, in combination with consideration of the above mentioned 
diagnostics. 

As mentioned in section 3.4, improper parameterization of <j) q and therefore 
iff could be possible sources of error in the parameterization of C E in stable conditions. 
Therefore, a modification of the (j) q and y/ q parameterizations will be investigated. 
Additionally, could the assumption of w q' w ~ 0 be incorrect? The results from section 

4.3 suggest yes. If the LI-7500 is only measuring wq , the inclusion of w q' w would 
increase the downward flux of total water and , perhaps, improve the comparison. 

In section 3.2, the roughness length, z , was introduced, yet no analysis of its 

effect on the parameterization of C E or (j) q was investigated in this case study. In light of 

the amount of work yet to be done to improve the overall quality of the analysis, we felt 
this was premature. However, once this work is completed, future studies will 
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incorporate how estimates of z q compare with TC 3.0 using improved estimates of 
<t> q and Y q • 

The proper parameterizations of the LHF and SHF by models such as CO AMPS 
are very important, especially in a stable boundary layer. These flux parameterizations, 
in addition to the surface stress parameterization, directly impact moisture, temperature 
and wind speed, near the surface. This in turn affects the higher level boundary layer 
structure including fog, clouds, and the low-level jet. In addition, fog development is 
very sensitive to SST which cools down and dries the low-level air under stable and 
saturation surface conditions. Since the stable surface layer is usually shallow, the errors 
in the LHF parameterization will greatly impact the saturation condition in the low levels. 

W 

From figure 1.1, it can be seen that TC 3.0 overestimates LHF by as much as 10—- in 

m 

stable conditions with fog. While the cumulative effect on the heat balance may be large, 

W 

it is unclear how much of an effect a 10—- error in the LHF would have on the SST over 

m 

shorter time scales. Fog will also likely be sensitive to the stress which generates 
turbulent mixing thereby effecting fog development (Wang 2004). 

In conclusion, further studies need to be conducted under stable conditions with 
fog to gain sufficient insight to develop improved parameterizations of scalar fluxes. 
Ideally, these investigations should be a combination of observational and modeling 
studies. Such an investigation has recently begun between WHOI and NRL Monterey. If 
an improved version of the TC 3.0 parameterization of the sensible and latent heat fluxes 
become utilized in COAMPS, there should be a significant improvement in the 
forecasting ability of the model. In turn, this will improve the Navy's ability for safe 
navigation and more accurate and effective mission planning. 
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APPENDIX A 


The discussion of the expansion of equation (3.4.14) from Edson (2004) is given 
by Businger (1982) and depends on the constraint of zero vertical mass flux for dry air: 

p a w = 0 

He uses this accepted governing constraint to derive an expression for the total mass flux 
(assuming no liquid water): 

pw = p a w+p v w= p(l-q)w+ p v w = p v w~ pw q 



pwqL e = pq (c pv - c pw ) wT + pL j^ 1 — 

7 7 

PwqCpJT = pqc pw wT' + c pw pT ^ "jV 

~pwKi = P~^=\ 

v~ c n 
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APPENDIX B 


This expansion of equation (3.4.32) from Edson (2004) begins by considering the 
mass flux. With liquid water included, the mass flux becomes: 

p w = p a w + p v w + p w w = p{l-q-q w )w+ pqw + pq w w 

As in Businger (1982), no contribution from horizontal convection is considered: 

p a w = p(l-q-q w )w=0 


pw= p v W + p w W = 


p 


( 1 - 0 - 0 *) 

The expansion of terms in (3.4.32) now become: 


(wq + wq ' w ) = p a w'r + p a w'r' w 


pw(\-q- q w ) c d T = pc, w'T'(\-q-q w ) 


pw{q + q w )c T = pc (q + q w )w'T' + pc j _ [w q + wq ' w ) 

l 1 -9-9,) 


pwqL e = pq (c -c )w'T'+ pL- --- -(w'q' + wq'J) 

P wb m 3 =Pp ---r(wV+wV')/> w 
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